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Abstract 

Recently, in an attempt to study disordered criticality in Quantum Hall sys- 
tems and d-wave superconductivity, it was found that two dimensional random 
Dirac fermion systems contain a line of critical points which is connected to 
the pure system. We use bosonization and current algebra to study properties 
of the critical line and calculate the exact scaling dimensions of all local op- 
erators. We find that the critical line contains an infinite number of relevant 

operators with negative scaling dimensions. 
71.10.-hx,71.28.-hd, 71.30.-Fh,11.40.-q 
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The transitions between quantum Hall (QH) states induced by disorder have been a long 
standing problem . Recently, it was pointed out that the transitions can also be induced by 
periodic potentials. The critical properties of these transitions in pure systems can then be 
studied through expansion and perturbative expansion 0,^. In particular, in absence 
of interactions, the pure critical point at the transition (between integer QH (IQH) states) 
can be described by an effective theory for Dirac fermions in 1+2 space-time dimensions. It 
is then natural to investigate what happens to the pure critical point if disorder is included. 
This point of view was taken recently by Ludwig et al. @]. 

In this letter we plan to study some interesting effects of impurities on the pure crit- 
ical point of IQH transition. Our starting point is the Dirac Hamiltonian in two spatial 
dimensions for the low energy spectrum of the pure system at criticality: 

Here, /i = 1,2 denote the two spatial directions and (71,72,75) = (^^1,(^2' ^3) the Pauli 
matrices. Three types of disorder are important: 

They correspond to random gauge potential, random chemical potential, and random mass. 
(The constant mass ^Qj^Mq is the parameter that controls the transition in the pure sys- 
tem 1^.) The strength of the impurities is parameterized by three positive parameters 
Qav.mi respectively, assuming the fluctuations of A^, V(x), and M(x) to be given by 
exp— I / (Px^Ajj^ + V^ + M"^). To leading order in the impurities strength, the three random- 
ness represent marginal perturbations to the pure critical point H^. Including the second 
order terms and in presence of only one type of impurities, it was shown in Ref. 0] that 
is marginally irrelevant, gy is marginally relevant, and gj^ is exactly marginal to all orders 
and generate a line of critical points starting at the free Dirac fermion model. It is the prop- 
erties of this critical line that we are going to concentrate on in this letter. The main result 
obtained here is that all the critical points on the critical line (except for the free Dirac 
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fermion model with g^^ = 0) contain infinitely many operators with negative dimensions. 
Those operators can be generated as the higher order terms in a non-Gaussian distribution 
of the random gauge potentials. Thus the critical points on the critical line are unstable in 
a very special sense - they have infinitely many relevant directions. 

The same Hamiltonian H = Hq + -f^imp ^"^^ i^s generalizations obtained by allowing for 
internal symmetries such as spin or isospin describe many non-interacting two-dimensional 
electronic systems characterized by isolated Fermi points in the presence of static disorder. 
Examples are degenerate semi-conductors, two-dimensional graphite sheets, tight-binding 
Hamiltonians in the flux phase, and dirty d-wave superconductors in two- dimensions 

Some exact properties of the critical line have been studied through bosonization in the 
replica approach In this letter we will use the supersymmetric formalism for non- 

interacting disordered systems 0. The advantage of the supersymmetric approach is that 
the operator content of the critical points can be easily obtained. This allows us to study a 
large class of local operators that may appear in the low energy effective theory. 

In the supersymmetric formalism, one begins by representing the unaveraged one-particle 
Green function at frequency by a path integral over pairs of Grassmann spinors ip, ip and 
pairs of complex spinors v^: 



Integration over (p and (p gives the denominator in the usual path-integral representation 
of Green functions. The infinitesimal positive number O"*" insures the convergence of the 
bosonic path integral. The impurity averaged (retarded) Green function is then given by 
the RHS of Eq. (|l|) with replaced by 




(1) 



where 



C^ = ^p i±)iuj-H)+iO+ + ^ {±){uj- H)+iO+ ip. 



[^7/.^^ + ^^ + ^ + (p [il^d^ + UJ + iO+J ip 
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+ ^ Y" (^^sV' + <^75¥')^ • 
Impurity averaging has turned a non-interacting problem into an interacting one. The 
effective action has an internal U{1/1) graded symmetry (supersymmetry) . 

The effective theory can be solved exactly if a; = gi^^ = Qj^ — gy — In this limit, the 
effective action defines a supersymmetric generalization of the Thirring model that we solve 
exactly in two ways. One approach is based on the current algebra with U{1/1) graded 
internal symmetry. It is quite general. We apply it to the particular case of the U{1/1) 
Thirring model. Another approach is specifically suited to the Thirring model, where it 
is possible to decouple the vector impurity from the spinors. In both cases the operator 
content on the critical hne is calculated. This will allow us to assess the stability of the hne 
of critical points. We begin with the latter approach. 

The partition function for the U{1/1) Thirring model is equivalent to averaging the 
one-particle Green function over impurities coupling to the currents only 



M 

X exp { - V57^ [d^ - + ^^$2)] ^ 

Here, we have rewritten the random gauge potential ^4^ as — + 9^$2 (^/i — ^^u^n^ 
€^2 — ~^2i — !)■ We have also reexponentiated the Jacobian Det 9^ for this change of 
variables with the help of two ghost fields a and a so as to preserve the unity of the partition 
function {A4 normalizes the measure of the impurity). One can decouple the spinors from 
the impurity potentials through a redefinition of the spinor fields. In terms of the chiral 
components ip± = |(1 ± Jb)^^, the decoupling transformation is: 



This transformation leaves the measure unchanged due to the super symmetry. We thus map 
the (interacting) U{1/1) Thirring model onto three independent and free sectors 



where we use complex coordinates z = Xi +ix2 and z = xi — ix2- Each sector is conformally 
invariant. The impurity strength does not appear explicitly, because it is hidden in the 
rotation from the original fields to the new ones. 

Any local operator in the original variables can be rewritten in terms of the four primary 
fields ^/'^ , ■ ■ ■ , and the two primary fields e**-*v^)*i and e*v^*2. Whereas the scaling 
dimension of e*^^*2 is positive, the scaling dimension of e*''*v^)*i is negative due to the 
non-compactness of the decoupling transformation. This in itself is not surprising. For 
example, primary fields with negative dimensions are common place in the treatment of 
strongly correlated electronic system [Q. They can turn a marginal interaction such as 
the Umklapp term into a marginally relevant interaction as the current-current interaction 
increases beyond a critical value 0. However, in the context of our supersymmetric model it 
leads for any given value of g^> to the existence of infinitely many local operators which 
are compatible with the U{1/1) symmetry and which carry negative dimensions. This is 
in sharp contrast to unitary theories (say the U(l) Thirring model) where operators with 
negative dimensions are not allowed. 

To see this, we begin by calculating the two point functions. The only non-vanishing 
ones are 




X e 





(2) 



as follows immediately from, say. 
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(4W±(0)) = (4(a:)^±(0)) 



(3) 



At the level of the two-point function, the scaling properties of the and the cp are not 
changed by the randomness of the gauge potential. However, this cannot not be the case for 
two-point functions of composite operators. For example, consider the composite operator 



O = ip^ ■■■ xy9^^ which is defined through point splitting. We find, with a 



calculation along the lines of Eq. (|]) , that O scales like ^0 

{0^{z,z)0 (0)) oc z-^''z-^~\ 

where the two conformal weights {h, h) are given by 

1 



(4) 



h = - 
2 
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m 



Here /i^2 are purely imaginary and real integer functions, respectively: 

/i = —i {rrii — m2 — + + — mg — -|- mg) , 

/2 = — {m^ + 7712 ~ "^3 ~ + ^5 + "^6 ~ ^7 ~ ■ 

Negative dimensional operators are now easily obtained. We consider the composite 
operator \E'„ „ defined by 



( n, n„ 


^1 


>o. 


712 


>0, 




^1 


<o. 


7l2 


>0, 






>o. 


7l2 


<0, 




^1 


<o. 


7l2 


<0, 



(5) 



and apply the result above. The conformal weights of these operators are 



u 1| I , fi'A T 1, I , S-A 

2' ^' 27r ^ 2' 2' ^' 27r ^ ^ 



(6) 



This demonstrates that for any given value of Qj^, there are infinitely many local composite 
operators with negative conformal weights. This is very different from the U{1) Thirring 
model which is a unitary theory and thus cannot support operators with negative dimensions. 
Furthermore, operators like are generated in the effective supersymmetric model if 

non-Gaussian moments in the probability distribution of the M{x) or V{x) are present. 

We now formulate a current algebra description of the model Eq. (|]) along the critical 
line oj = Qj^j = (7y = 0. This provides us with an independent check of Eq. (|), and more 
importantly, can be generalized to any conformally invariant theory with internal U{1/1) 
super symmetry. We first construct the current algebra description for the free model, and 
then show that the current algebra approach also applies to any positive value of Qj^. 

When = we have the free theory 

It describes two decoupled sectors labeled by the ± subscripts of the chiral components. 
One can study this theory by canonically quantizing it along contours of equal "time" in 
the x^-X2 Euclidean plane. In what follows, we shall use, without loss of generality, equal 
X2 lines. The canonical momenta associated with the fields il)± are 

TT ^-^0 ,t TT t 

= , N = =Fm, n = — — — - = =Fv?±, 

d[d2il)±) o{d2(p±) 
and the equal ("time") quantization conditions are {ip±{xi,X2), 11^^ (x'^, X2)} = i6{xi — x[) 

and [(p±{xi, X2), Il^^{x[, X2)] = i6{xi —x[). The time-ordered Green's functions for ip, ip can 
be obtained from the classical equations of motion found from the Lagrangian, together with 
the above commutation relations used as boundary conditions. One recovers the two-point 
functions of Eqs. (0) obtained from the path integral. 

The free model has an internal graded symmetry f/(l/l), generated by two bosonic and 
two fermionic currents. The bosonic U{1) symmetries are {t/j^,(p^) (e^^^ip^^e^^^fp^). 
The fermionic symmetries are generated by the infinitesimal transformations {6ijj^,6(p^) = 
i6<^±, —0'ip±), where 6 and 6 are Grassmann variables. The normal ordered conserved cur- 
rents associated with these symmetries are: 
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jt{xi,X2) = ip±{xi + e,X2)ilJ±{xi - e,X2) - 

J^{Xi, X2) = (f±{xi + €, X2)(f±{Xi - e, X2) + 

r]±{xi,X2) = <-P±{xi + e,X2)iJ±{xi - e,X2), 
fj±{xi,X2) = ^±{xi + e, X2)(p±{xi - e, X2), 

where point splitting has been used (e ^ 0). It is convenient to define the currents J± — 
J% + J± and j± — J% — J'^. Using the equal X2 commutation relations for the ip,(fi^s, we find 
that the non-vanishing commutators are: 

[J±{Xi,X2),j±{x[,X2)] = ±-5'{xi -X[), 
[j±{XuX2),V±{x[,X2)] = ±2i S{Xi -X[) V±{x'i:X2), 

[j±{xuX2),fi±{x[,X2)] = T2i S{xi - x[) fi±{x[,X2), 
{V±{Xl,X2),fj±{x[,X2)} = Ti[^S'{xi -x[) 

+ S{xi- x[) J±{x[,X2) 

All commutators between the + and — sectors vanish. The chiral components and J_ 
are the holomorphic and antiholomorphic components of Ji = J+ + J_ and J2 — i{J+ — J-), 
respectively. 

From the equations of motion derived using Cq, we find that fields in the + (— ) sectors 
such as ipj., (/?_!_, and = J^,j^,r]^,fi^, depend only on z {z). Such fields are called 
chiral fields. The same result can also be obtained from the quantum equation of motion 
820 = i[H^ 1 O], where is derived from our canonical quantization. However, it is more 
convenient to express as a bilinear form in the currents. We find that = + 
with 

Hl = iTi j dxi {J±J± + J±j± + r]±fj± - f}±ri±) (7) 

exactly reproduces the equations of motion for the currents and ■0, V- It then follows that, 
for chiral fields, say in the + sector, Oj^{z) = e'^^+^O_^_{0)e~^+^ . Notice that is invariant 
under U{1/1) transformations. 



The impurity averaging induces an interaction C\-at = i^-Jf^ = '^2(7^J+J_ in the La- 
grangian. This interaction, in turn, corresponds to TYint = —i'^QA J+J-- The new Hamilto- 
nian can be diagonahzed, while keeping the + and — sectors independent, by a redefinition of 
the currents J± = J±, 'U± = rj±, fj^ = fj±, j± =j± — ^J± — ^J^. For convenience, we have 
chosen to maintain the j charge density unchanged, i.e., j2 = — j-) = — 3-) = 32- 
The chiral components oi H = J dxiH are now given by 

= ( 1 + ^ ) J± J± + Jj± + ^±^± - ^±77±. (8) 

With this redefinitions of the currents, one can check that the only current commutator 
that changes is \3±{x\,X2),3±{x'y,X2)\ = Ti-^ 6'{xi —x\). We thus see that the interacting 
theory also decouples into a holomorphic and antiholomorphic sector. The Hamiltonian 
Eq. (P) and the commutators between 7!^ also allow us to calculate the exact A^-point 
correlations between the currents. One can show that all the correlations have algebraic 
decay, which implies that the interacting model describes a critical point. This demonstrates 
the conformal invariance of the theory for arbitrary values of g^. 

We now sketch how the current algebra approach can be used to calculate the anomalous 
dimension of a generic operator O constructed locally from powers of ip^s and (p. Due to the 
conformal symmetry of the critical point, the two-point correlation function on the plane, 
Eq. @, can be mapped onto the one on a cylinder: {O{z,z)O{0)) ~ ^thz-zhz large X2- 
We have taken the circular direction of the cylinder to be < Xi < 27r. From 

{O{z,z)O{0)) = (0|C(0)e-^^++^"^-O(0)|0), (9) 

we see that the large X2 behavior is controlled by the state \0). Here, we have expanded 
C(0)|0) into the eigenstates of iH [{iHY = iH], and \0) is the one with the minimal 
eigenvalue Eq^ + Eq~ , where H±\0) = Eq^\0). The dimensions of O are then simply 
given by {h,h) = {Eq+,Eq-). 

Notice now that the Fourier components of the currents, 1^.^ = / d^i I±{xi)e'^"'^^ , satisfy 
[iH, I±-n] = ='=^-^±;n- Heucc, the modes I±.n raise and lower the eigenvalues of iH. One can 

9 



then show that the state \0) is annihilated by the lowering operators: /|._„|(9) = 
0, for n > 0, so that using Eq. (H), 
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2 



9 a- 



lH±\0) = - il + ^)J±;oJ±;0 + J±;Oj± 
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+ ^±;0^±;0-^±;0^±;o]|C)- (10) 
By reexpressing Eq. (|r^) in terms of the unrotated currents, one finally obtains 

h\0) = /ifreelO) - f^J+;oJ-;o|0), (11) 

where hf^ee is the weight of O in the absence of interactions. A similar expression can be 
found for h. Thus we reach the very simple result that the conformal weights {h, h) of O 
can be calculated from the charges of O for unrotated currents. 

The J charges of operators constructed from the ip, y^'s can be easily obtained from 
their commutators with J±-fl. The conformal weights for ^n^n^ Eq. then follows 
immediately by recognizing that its unrotated charges are (J+;o, <^-;o) = {.—ni,n2) and that 
its free conformal weights are (/ifree, ^free) = (^l^il, ||^2|)- 

We believe that negative dimensional operators are generic to disordered critical points, 
as one can infer from this work and the 2 + e expansion |TT|. In general, negative dimensional 



operators are allowed since supersymmetric models for the disordered critical points are not 
unitary and, as such, contain states with negative norm. Therefore, in the search for physical 
disordered critical points, it is crucial to find under what conditions the generally allowed 
negative dimensional operators do not spoil the critical points. 

We would also like to point out that the supersymmetric model studied here can be 
viewed as the fermionized disordered X—Y model studied in [|12||. Our results are then very 
suggestive of a similar instability that might occur along the critical line of the disordered 



X — Y model [0. A more detailed study will appear elsewhere. 
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